Elastic nonlinearities in a one-dimensional model of fracture 
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The dynamics of rapid brittle cracks is commonly studied in the framework of linear elastic fracture 
mechanics where nonlinearities are neglected. However, recent experimental and theoretical work 
demonstrated explicitly the importance of elastic nonlinearities in fracture dynamics. We study 
two simple one-dimensional models of fracture in order to gain insights about the role of elastic 
nonlinearities and the implications of their exclusion in the common linear elastic approximation. 
In one model we consider the decohesion of a nonlinear elastic membrane from a substrate. In a 
second model we follow the philosophy of linear elastic fracture mechanics and study a linearized 
version of the nonlinear model. By analyzing the steady state solutions, the velocity-load relations 
and the response to perturbations of the two models we show that the linear approximation fails at 
finite crack tip velocities. We highlight certain features of the breakdown of the linear theory and 
discuss possible implications of our results to higher dimensional systems. 



I. INTRODUCTION 

The dynamics of rapid brittle cracks exhibits a rich 
phenomenology that is not yet well understood. For ex- 
ample, crack tip instabilities (such as a side-branching 
instability [l| and an oscillatory one Q), that were 
shown to govern fracture dynamics at high propagation 
velocities, are poorly understood from a fundamental 
point of view. The major stumbling block in devel- 
oping a well-established theory of these phenomena is 
our lack of understanding of the physics of the "frac- 
ture process zone" within which nonlinear deformation, 
dissipation and material separation processes take place. 
The classic approach of linear elastic fracture mechanics 
(LEFM) assumes infinitesimal deformation outside this 
typically small process zone and predicts asymptotically 
"diverging- like" strain and stress fields [3j. Under these 
assumptions the energy flux into the process zone is cal- 
culated and an equation for the rate of crack growth is ob- 
tained by equating it to an unknown dissipation function 
r. This quantity lumps together all the poorly under- 
stood nonlinear and dissipative properties of the process 
zone dynamics Q. 

Very recent experimental and theoretical work demon- 
strated explicitly the existence of a nonlinear elastic zone 
in the near vicinity of a rapidly moving crack tip, where 
the deformation fields were shown to be quite different 
from those predicted by LEFM 0, Q • One implication 
of these findings is that although LEFM may provide 
reasonable estimates of the energy flowing to the pro- 
cess zone, it fails to represent properly the ways in which 
breaking stresses are being transmitted to the crack tip. 
Therefore, as long as the path of the crack is known and is 
stable against perturbations, the energy based approach 
of LEFM seems useful; however, as such conditions are 
rarely met and the question of path stability is usually 
of prime importance, the near tip deformation and stress 
fields themselves may play a central role in describing 
fast fracture and the associated instabilities. Moreover, 
these tip instabilities seem to involve a non-geometrical 



lengthscale (for example, the minimal side-branch length 
[f||3] or the wavelength of oscillations 0,11]) that is miss- 
ing in LEFM. Thus, the findings of 0, that suggest 
the existence of a dynamical lengthscale associated with 
the nonlinear elastic zone, shed new light on the search 
for a missing lengthscale. 

Motivated by these recent results, we aim at gain- 
ing additional insights about the possible roles played 
by elastic nonlinearities in fracture dynamics and about 
the possible implications of neglecting these effects in the 
common linear elastic approach. For that aim we study 
in this paper simple one-dimensional models of fracture 
in which an elastic membrane is being detached from 
a substrate by the propagation of a decohesion front. 
The models are being defined by the following differential 
equation for the scalar deformation u(x, t) of the mem- 
brane 



pdttu = d x s - k 2 (u- 5) - 4>(u,d t u) 



(1) 



where p is the linear mass density of the membrane. The 
first term on the right-hand-side represents the force per 
unit length due to the deformation u(x,t), where the 
stress s(x,t) is related to u(x,t) through a constitutive 
law of the form 



s = pT(d x u) 



(2) 



Here p is the elastic modulus of the membrane and the 
functional T{d x u) represents a general nonlinear stress- 
strain relation. The strain is the displacement gradient 
d x u. The second term on the right-hand-side represents 
the loading of the membrane by elastic springs whose 
spring constant is k 2 and whose natural length is reached 
when u(x, t) = S. Note that p has the dimension of force, 
while n 2 has the dimension of force/squared length. The 
term <f)(u, dtu) represents the visco-elastic interaction of 
the membrane with the substrate. This term should also 
include a criterion for detachment from the substrate, 
serving as a fracture criterion. 

Similar one-dimensional models were studied previ- 
ously in relation to various aspects of fracture dynamics, 
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see for example 0, [13, HH, EH ■ Our model follows the 
spirit of these works and its linear app roximation is simi- 
lar to the model discussed in Ref . [Ty] . Our strategy is to 
study two related models with the same (f>(u, dtu), where 
in one we consider the solution for a nonlinear T{d x u) 
and in the other we follow the philosophy of LEFM and 
linearize J-{d x u) first and then solve the model. Our 
main goal is to compare various aspects of the dynamics 
of the nonlinear and linearized models in order to gain 
insights about the role of elastic nonlinearities in the dy- 
namics of fast brittle cracks. 

Our results show that for sufficiently large k there ex- 
ists a range of small crack tip velocities v (i.e. «Cc, 
where c= \J [ij p is the linear wave speed), such that the 
linear theory provides reasonable approximations to the 
nonlinear theory. However, as the velocity v increases 
the linear approximation deteriorates progressively until 
it fails to capture important aspects of the dynamics. For 
smaller values of k there exists no range of validity for the 
linear approximation. Our conclusion is that in our sim- 
ple one-dimensional model of fracture LEFM inevitably 
breaks down at finite crack tip velocities, in agreement 
with the findings of 0, [f| . This breakdown can be of sig- 
nificance to the understanding of crack tip instabilities in 
higher dimensional and more realistic fracture problems. 
Specifically, we show how elastic nonlinearities affect the 
limiting crack velocity, the strain and stress fields near 
the moving crack tip and the dynamical time and length 
scales involved in the physics of the near tip region. 

In Sect. [IT] we present the models by first introducing a 
nonlinear constitutive law that can capture both soften- 
ing and stiffening behaviors. Then we present a linearized 
version of this nonlinear model. In Sect. IIIII we solve for 
the steady states of the models and discuss the result- 
ing deformation profiles, the load-velocity curves and the 
limiting crack tip velocities. In Sect. IIVI we study the re- 
sponse of the steady states to perturbations in the crack 
tip location within the framework of a linear stability 
analysis. In Sect. [V]we summarize the results and dis- 
cuss their possible implications to more realistic, higher 
dimensional, systems. 



II. THE MODELS 

In order to complete the definition of the models to 
be considered below, we should supplement Eq. (JTJ) with 
a constitutive law T(d x u) in Eq. ([2]) and to specify the 
form of the interaction of the membrane with the sub- 
strate 4>{u,d t u). For the latter we choose 



(u, dtu) — a 2 uH(uo — u) + rjdtuH(uo — u) 



(3) 



This is the cohesion force that binds the membrane to 
the substrate and can be thought of as the action of 
visco-elastic springs with a spring constant a 2 , a small 
friction/ viscosity-like coefficient rj and a breaking thresh- 
old uq. The Heaviside step function H(-) represents the 
irreversible breaking of the cohesion springs when the 
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FIG. 1: (Color online) A sketch of the one-dimensional model. 
A membrane, whose profile is given by the deformation u(x, t), 
propagates from right to left (in this example) at a velocity 
v. The vertical solid lines represent the loading springs. The 
vertical dashed lines represent the cohesive springs, uo is the 
breaking threshold of the cohesive springs, <5 is the deforma- 
tion of the membrane for which the loading springs are com- 
pletely relaxed (see Eq. JTJ) and u_ is the asymptotic value 
of the deformation in the negative x direction (see Eq. (|14p). 



displacement exceeds uo- The small dissipation associ- 
ated with r] is essential to ensure the existence of steady 
states, see [1(3] and below for details. A sketch of the 
one-dimensional model is shown in Fig. [TJ 

The constitutive law for the membrane is chosen to be 



(1 + d x u) - (1 + d x u)~ 



(4) 



This law corresponds to the uni-axial behavior of the 
tensorial neo-Hookean constitutive law considered in [f| . 
For small strains, Eq. ((4]) can be linearized, yielding 



F{d x u) ~d x u + [(d x u) 2 ] for d x u < 1 



(5) 



The constitutive law of Eq. J3J) is plotted in Fig. [2] For 
positive strains, d x u > 0, the constitutive law exhibits 
nonlinear softening, i.e. the local tangent to the curve 
is smaller than the tangent at infinitesimal strains. For 
negative strains, d x u < 0, the constitutive law exhibits 
nonlinear stiffening, i.e. the local tangent to the curve is 
larger than the tangent at infinitesimal d x u. Note that 
in the nonlinear softening case the tangent approaches 
1/3 at large strains. Usually, fracture is associated with 
positive strains. However, in our simple model the sign 
of d x u is determined by the direction of crack propaga- 
tion; when the crack propagates from right to left we have 
d x u > and the softening branch is selected, while when 
it propagates from left to right we have d x u < and the 
stiffening branch is selected. Thus, both nonlinear be- 
haviors can be incorporated into our model and we refer 
to these cases as the nonlinear softening and nonlinear 
stiffening models respectively. When the linearized rela- 
tion of Eq. is used instead of Eq. ^ , we refer to the 
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for the nonlinear models and 




FIG. 2: (Color online) !F = s/fi as a function of d x u (solid 
line). The linear approximation at small d x u is added (dashed 
line). The branch corresponding to d x u>0 exhibits nonlinear 
softening, while the one that corresponds to d x u<0 exhibits 
nonlinear stiffening. 



model as the linear one. Note that in that case the direc- 
tion of crack propagation is irrelevant as the linearized 
relation is symmetric. 

In order to prepare Eq. |T]) for the analysis to come 
we nondimcnsionalize all of the quantities by measuring 
length in units of uq, velocity in units of c and force in 
units of fx. For simplicity we denote all the nondimen- 
sionalized quantities using their original notations. Fur- 
thermore, we set a = l. We thus end with the following 
equation 

duu — d x T{d x u) — K 2 (u — 5) — uH(l—u) — ridtuH(l — u) , 



where 



d x F(d x u) 



(i + d xU y 



(6) 



(7) 



d x T(d x u) ~ d xx u 



(8) 



for the linear one. 



III. STEADY STATE SOLUTIONS 

Our goal in this section is to solve for the steady states 
of both the linear and nonlinear models and to compare 
various properties of these solutions. The tip of the crack 
is defined as the point where the cohesive springs break 
and assumed to propagate at a constant velocity In a 
co-moving frame x' = x + vt the crack tip lies at the origin 
and u(0) = 1. We now look for a steady state solution of 
the form u(x + vt). We start by considering the linear 
model, for which one can obtain analytic results. 

A. The linear model 

For steady state conditions, Eq. flB}, with Eq. ©, 
reduces to an ordinary differential equation of the form 

= (1 - v 2 )d xx u - k 2 (u -6)-(u+ vr)d x u)H{-x) , (9) 

where we replaced x' with x for the simplicity of the 
notation. The appearance of a step function, modeling 
the rupture of the cohesive springs at a critical displace- 
ment, implies that the differential equation is ill-defined 
at x = 0. Therefore, we treat separately the positive and 
negative x domains and demand continuity of u(0) and 
d x u(0). We first consider the domain x>0. In the limit 
x^+oo we have u— >6. Therefore, we assume a solution 
of the form u = 6 +(1 — 6) exp (kx), that satisfies u(Q) = 1. 
Substituting in Eq. @, we obtain a simple second order 
algebraic equation for k. We choose the negative root 
since x>0 (the solution must be bounded). Following a 
similar procedure for x < 0, we arrive at 



u(x) = 6 + (1 — 6) exp — 



u(x) 



5k 2 



1 + K 2 



1 - 



6k 2 



1 + K 2 



exp 



rjv + y/rjhP + 4(1 + n 2 )(l - v 2 ) 
2(1 -v 2 ) 



for x > , 
, for x < . 



(10) 



where the velocity v is still undetermined. The velocity 
is determined by demanding that c^i^O) is continuous. 
Note that higher order derivatives at x = are discontin- 
uous due to the discontinuity of the force at this point. 



Using the continuity of the derivative at x = we obtain 



V(S,K,T]) 



6 2 - 6 2 



^/(6 2 - 6 2 ) 2 + k 2 t] 2 (8 - 1) 2 (6 - 6 2 ) : 



, (11) 
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with 



S c = VI 



(12) 



We are now able to interpret the steady state solutions 
of the linear model in terms of standard LEFM concepts. 
<5, with a fixed k, can represent the load and r\ is a mate- 
rial parameter that is related to dissipation. Within this 
interpretation, Eq. (jTTJ) tells us that S c is the critical load 
needed to initiate crack propagation. It is a direct ana- 
log of the Griffith criterion 0. Furthermore, Eq. (TTTj) 
with a fixed K represents the velocity-load relation of the 
model. For 5 close to 5 C , the velocity of the crack is 
much smaller than the wave speed, while for sufficiently 
large 6 the tip velocity approaches the wave speed. Thus, 
the limiting crack velocity is the wave speed, in complete 
analogy with the common prediction of LEFM in higher 
dimensions, where the limiting velocity is the Rayleigh 
wave speed 

The role of the small dissipation coefficient 77 is most 
clearly demonstrated by multiplying Eq. ^ by d x u and 
integrating from —00 to 00 to obtain 



1 



(1 



rjv 



1 



(d x u) 2 dx = -k\5-u-Y , (13) 



where 



5k 2 



1 + K 2 



(14) 



is the value of«asn-oo. It is important to note that 
this result holds for any elastic constitutive law T(d x u) 
and not only for a linear one. The first term on the 
left-hand-side of Eq. (fT3|) is simply the energy per unit 
length needed to break the cohesive bonds. Therefore, it 
is simply the bare surface energy 7 



1, 9 X 

7= 2(1-"-) • 



(15) 



The second term on the left-hand-side is the fric- 
tional/ viscous dissipation. Together, these two terms re- 
sult in the fracture energy T(v) [3j 



T(v, 77) = 7 + rjv I (d x u) 2 dx 



(16) 



The right-hand-side of Eq. (Ti~3"|) is the energy per unit 
length stored far ahead of the crack tip. Therefore, this 
is the so-called "energy release rate" G(v) [|[ 



Giv) = r {d - u - )= wT^) 2 



(17) 



where S(v) is obtained by inverting Eq. (TTTj) . Therefore, 
we can rewrite the energy balance of Eq. ([13]) in the com- 
mon LEFM form G(v) = T(v,T]) Q. By reexamining Eq. 
(Till)) we observe that the velocity v is coupled to 77 such 
that when 77 vanishes v does not appear in the equation. 



Therefore, for 77 = 0, there exist steady states only when 
S = S c and in this case, since G{v) — T(v, 0) indepen- 
dently of v, the crack can propagate at any velocity. For 
5>S C we have G>T and there exist no steady states. In 
that case the crack is expected to accelerate toward the 
limiting velocity v = 1 . This result is consistent with the 
77 ^0 limit of Eq. (TTTj) . However, one should be cautious 
as this limit predicts that v = 1 also for 6 < S c , which in 
light of the discussion above, is wrong. This observation 
was made previously for a related model [lof . 

In order to set the stage for the comparison between 
the linear and nonlinear models below, we should ask 
what one can learn about this comparison from the 
steady states of the linear model alone. Intuitively, it 
is quite clear what determines the range of validity of 
the linear approximation: it is expected to breakdown 
if strains significantly larger than a few percent develop 
over a large enough region near the crack tip. Mathe- 
matically speaking, we expect the linear approximation 
to hold if (d x u) 2 is sufficiently smaller than d x u almost 
everywhere. In order to make this qualitative observa- 
tion more quantitative, we note that the largest strain 
derived from Eq. (JTDJ) , which occurs at a; = 0, is 



d x u(0) 



(5-1)k 



(18) 



This quantity can be interpreted as the analog of the 
stress intensity factor of LEFM [3j, since it provides a 
measure of the typical strains near the tip of the crack. 
Moreover, it is an increasing function of the load 5 and 
k, and more importantly, it diverges in the limit v — > 1. 
Thus, we already learn that even if there exists a low ve- 
locities range where the linear approximation holds, then 
it will breakdown at a finite, possibly high, velocity. In 
fact, the issue of whether there exists a low velocities 
range of validity of the linear approximation can be fur- 
ther elucidated by substituting 6 — 5 C and v = in Eq. 
(fig)) to obtain the following inequality 



d x u(x;v = 0) < (y/l + k- 2 - 1)« < d x u(x = 0;v) . (19) 



This inequality suggests that the value of (vl-j- kt 2 — 
1)k, in comparison to unity, determines whether there 
exists a range of low velocities where the linear approx- 
imation is valid or not. For example, for k — 2 we ob- 
tain (Vi + k 2 — 1)k~0.24 and we expect the linearized 
model to provide a reasonable approximation to the non- 
linear model at low velocities (see Fig. [2]), while for 
k = 0.3 we obtain (vT+ k~ 2 — 1)k ~ 0.74, for which 
we expect that no range of validity for the linear approx- 
imation exists. In the next subsection we solve for the 
steady states of the nonlinear models and present a de- 
tailed comparison with the results obtained above for the 
linear model. 



5 



0.8 



0.6 



0.4 



0.2 











* 

■* / 

/ / 
I / 
















// / 0.5 










* / / OA 










ft * > 0.3 










_ if* 0.2 










'J* 01 










¥ 1.0001 


1.0003 g^g 1.0005 

c 


1 .0007 









1.05 



1.1 



1.15 

8/8 

c 



1.2 



1.25 



1.3 



FIG. 3: (Color online) Main panel: The velocity v vs. the nor- 
malized load <5/<5 c for the linear model (solid line) - Eq. (Illfl . 
the nonlinear softening model (dashed line) and the nonlin- 
ear stiffening model (dashed-dotted line), all with k = 0.3 and 
r; = 0.3. The result for the linear model is given in Eq. 
Note that there exist no smooth steady states with v > 1 in 
the nonlinear stiffening model. In this range shock behavior 
is expected (not shown here). Inset: The velocity v vs. the 
normalized load 5/S c for the linear model (solid line) and the 
nonlinear softening model (dashed line) with k = 2 and r\ = 0.1. 



B. The nonlinear models 

Our aim in this subsection is to solve for the steady 
states of the nonlinear models and compare the results 
with those of the linear model. Substituting Eq. ([7]) 
into Eq. using the steady state assumption u(x±vt) 
and replacing x' = x ±vt with x for the simplicity of the 
notation, we obtain 



'1 2 
v 



d xx u + 



2 *i" *tjb. 

\3 ~ 



3(1 + cU) 3 
k 2 (u — 5) + (u =p vrjd x u) H (tx) ■ 



(20) 



We note that the =F signs correspond to nonlinear soft- 
ening and nonlinear stiffening respectively For the for- 
mer the crack propagates in the negative x-direction and 
d x u > such that the softening branch of the consti- 
tutive law of Eq. (j4]) is selected (see Fig. [2]), while for 
the latter the crack propagates in the positive x-direction 
and d x u < such that the stiffening branch of the con- 
stitutive law of Eq. ([4]) is selected (see Fig. [2]). The 
appearance of nonlinearities in Eq. (|20p entails a numer- 
ical solution. We solve the problem by using the shooting 
method. We first guess d x u(0) and v, and then integrate 
Eq. (|20[) from x — in both directions using a forth order 
Runge-Kutta integration scheme, where u(0) = l is used. 
We then improve the guess until u(x — > —00) — > S and 
u{x— >oo)— >(5k 2 /(1 + k 2 ) are approached monotonically. 

We first consider the velocity-load relations for both 
the linear and nonlinear models. In the main panel of 
Fig. [3] the propagation velocity v is plotted as a function 
of S/S c for the linear and nonlinear models with k = 0.3 



and r\ = 0.3. We first discuss the effect of elastic non- 
linearities on the limiting crack velocity. In the linear 
case (solid line), as discussed in relation to Eq. (TlT|) . 
the limiting velocity is v— * 1, which is the linear elastic 
wave speed. In the nonlinear softening case (dashed line) 
the limiting velocity is substantially smaller, 1/V3- 
To understand this, recall that the limiting crack veloc- 
ity is determined by the speed of small amplitude waves 
traveling near the tip of the crack, since these waves de- 
termine the rate at which energy is being transferred to 
the crack tip for breaking cohesive bonds. The speed of 
these waves, in higher dimensional models, is determined 
by the properties of the bulk material. In the simple 
one-dimensional models considered here there is no clear 
separation between the "bulk", that is represented by 
the elastic membrane, the external loading that is rep- 
resented by the n 2 (u — S) term and the cohesive force 
4>(u,dtu); therefore, the speed of small amplitude waves 
is affected by the loading and interfacial cohesion, in ad- 
dition to the "bulk" properties. However, for the sake of 
obtaining a physical understanding of the effect of elastic 
nonlinearities on the limiting crack velocity we consider 
in the discussion below only the elastic properties of the 
membrane. 

In the linear model the small amplitude wave speed 
is independent of deformation and equals to the small 
strains wave speed, therefore v — ► 1 in our dimension- 
less units. However, in the nonlinear models the small 
amplitude wave speed, that is determined by the local 
tangent to the stress-strain curve, depends on the state 
of deformation near the moving crack tip. In the soft- 
ening branch of the stress-strain curve presented in Fig. 
[3J the local tangent decreases with increasing strain un- 
til it approaches 1/3. Since the crack tip concentrates 
large strains, and the magnitude of these strains increases 
with increasing propagation velocities, the tip velocity 
is determined by the square root of the limiting local 
tangent 1/V3- In the stiffening branch of the stress- 
strain curve presented in Fig. [21 the local tangent con- 
tinuously increases with increasing the magnitude of the 
strain (the strain itself is negative in this case). There- 
fore, we expect cracks to propagate faster in the nonlinear 
stiffening model, compared to the linear model. More- 
over, if the crack propagates at velocities higher than the 
small strains wave speed, we expect the development of 
shocks. The velocity-load relation for the nonlinear stiff- 
ening model is shown in Fig. [3] (dashed-dotted line). As 
expected, for a given load, the propagation velocity is 
higher than in the linear and nonlinear softening models. 
Moreover, when the velocity approaches the small strains 
wave speed v— ► 1, we failed to find smooth steady states. 
This point marks the onset of shock development. Note, 
however, that we expect the existence of steady states 
with v > 1, thus the limiting velocity is expected to be 
higher than in the linear case. Similar ideas about the ef- 
fect of elastic nonlinearities on the limiting crack velocity 
were discussed previously in the literature [H, [13, . 

As discussed at the end of the previous subsection, we 
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FIG. 4: (Color online) (a) Main panel: The steady state 
strain d x u(x) for the linear (solid line) and nonlinear soften- 
ing (dashed line) models with 5=1.25<5 C , K = 0.3 and r/ — 0.3. 
Note that the velocities of the two profiles are different, see 
Fig. [3] Inset: The displacement u(x) in the near crack tip 
region. The arrow indicates the direction of propagation, (b) 
The same as (a), but for the linear (solid line) and nonlinear 
stiffening (dashed line) models with S — 1.058 c . (c) The same 
as (a), but with S — 1.05<5 C . (d) The same as (a), but with 
<S=1.000W C , k = 2 and 77 = 0.1 



expect the linear approximation to fail at all velocities for 
k = 0.3. This is demonstrated in Fig. [3]by the large differ- 
ence in propagation velocities. In fact, even in the limit 
6/8 c —> 1 the differences in the velocity of propagation are 
large due to the significantly different initial slopes of the 
velocity-load curves. In the inset of Fig. [3] we compare 
the velocity-load relations for the linear and nonlinear 
softening models with k = 2 and 77 = 0.1. For this value 
of k we expect the linear model to provide a reasonable 
approximation to the nonlinear model at small velocities. 
Indeed, the inset of Fig. [3] shows that at least as far as 
the velocity of propagation is concerned the velocities in 
the low v regime are similar in the two models, though 
they separate progressively with increasing load. 

We now turn to compare the deformation in the var- 
ious models. In the main panel of Fig. 0^ the strain 
distributions for the linear and nonlinear softening mod- 
els with 8 = 1.258c, k — 0.3 and 77 = 0.3 are shown. In 
light of the results shown in Fig. [31 we expect large dif- 
ferences in the strain near the crack tip for such a value 
of the load. Indeed, rather dramatic differences are ob- 
served near the tip of the crack, where the strain in the 
nonlinear softening case is significantly higher than in the 
linear case. This effect is in qualitative agreement with 
the findings of [3, IE]. In the inset, the displacement in 
the near tip region is shown. We note that according to 
Eq. (fT3")l the difference in the integral of (d x u) 2 over the 
negative x-axis determines the difference in propagation 
velocity when 8, k and 77 are fixed. The large difference 
in the near tip strains observed in Fig. [4^, is consistent 
with the large difference in velocity for <5 = 1.25# c in Fig. 
[3l In Fig. |4Jd the strain and displacement fields for the 
linear and nonlinear stiffening models with 8 = 1.05<5 C , 
k = 0.3 and 77 = 0.3 are shown. In this case the nonlinear 
strains are smaller than the linear ones. The correspond- 



ing comparison (i.e. for 8 = 1 .05<5 C ) with the nonlinear 
softening model is shown in Figs. 2)d and Hfc, both 
with 8 — 1.05(5 C exhibit smaller differences between the 
linear and nonlinear models compared to Fig. [^i where 
8 = 1.25<5 C , demonstrating quantitatively how the linear 
approximation deteriorates with increasing crack veloc- 
ity. Note, however, that even for 8 = 1.05<5 C the differences 
are non-negligible, in agreement with Fig. [3] that implies 
that there is no range of validity for the linear approxima- 
tion for k = 0.3. In Fig. |4ji the strain and displacement 
fields for the linear and nonlinear softening models with 
8 = 1.00018 C , k~2 and 77 = 0.1 are shown, cf. the inset of 
Fig. [3l For this larger value of k, the linear model pro- 
vides a reasonable approximation for the nonlinear one, 
where moderate differences in strain are observed only 
very near vicinity of the crack tip. 

In Fig. [5^ the elastic strain energy distributions for 
the three models with 8 = 1.05<5 C , k = 0.3 and 77 = 0.3 
are shown. The strain energy functional corresponding 
to the stress-strain relation of Eq. ^ is 



U = 



1 



l(l + d x u) 2 + j ' 



d x u 



(21) 



where Eq. ((¥]) is recovered using the following functional 
derivative 



8U 
8{d x u) 



(22) 



Note that in the common linear clastic approximation, 
Eq. (J3I]) yields 



U 



\{d x uf 



- O [{d x uf 



(23) 



In Fig. [SJd the kinetic energy distributions T = 
v 2 (d x u) 2 /2 for the three models with 5 = 1.05<5 C , k = 0.3 
and 77 = 0.3 are shown. The corresponding stress distri- 
butions are shown in Fig. [SJ3. Note that the magnitude 
of the stress at the tip is similar for the linear and non- 
linear softening models, but is different for the nonlinear 
stiffening case. This can be understood as follows: for 
the linear model we have s~d x u. For the nonlinear soft- 
ening model Eq. ((4]) is approximated as s ~ |(1 + d x u) 
for the large strains near the tip. The difference between 
these expressions is compensated by the larger strains in 
the nonlinear softening case, cf. Fig. yielding similar 
values for s. However, for the nonlinear stiffening model, 
Eq. ([J]) is approximated as s~ — ~(1 + d x u)~ 2 for the 
large negative strains near the tip, cf. Fig. [4Jd, resulting 
in a negative and somewhat larger in magnitude stress in 
this case. In spite of the fact that the details of the field 
distributions appearing in Figs. EMSt are specific to the 
simple models considered here, they give one a sense of 
the type of errors expected in models that exclude elastic 
nonlinearities even at moderate crack propagation veloc- 
ities. Possible implications of these differences for ques- 
tions of crack tip stability will be discussed in Sect. IIVI 
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FIG. 5: (Color online) (a) The steady state potential energy 
distribution U for the linear (solid line), nonlinear soften- 
ing (dashed line) and nonlinear stiffening (dashed-dotted line) 
models with <5 = 1.05<5 C , k = 0.3 and 7/ = 0.3. (b) The same as 
(a), but for the steady state kinetic energy distribution T. (c) 
The same as (a), but for the steady state stress distribution s. 
(d) The same as (a), but for the small amplitude local wave 
speed distribution c; oc of Eq. (|24[). In all panels the data for 
the nonlinear stiffening model was transformed according to 
x^> —x for the sake of comparison with the other models. 



The last issue to be discussed in relation to the steady 
state solutions concerns material lengthscales. It is well 
known that linear elasticity contains no intrinsic length- 
scale, while the instabilities of dynamics fracture indi- 
cate that some non-geometrical lengthscale is involved 
@) B 0j 111- The existence of elastic nonlinearities nat- 
urally suggests a lengthscale 0, Q , which is simply the 
size of the nonlinear zone. More precisely, a lengthscale 
can be defined as the size of the region in which material 
properties become deformation dependent, for example 
the region where the small amplitude local wave speed 
differs from the linear elastic wave speed. Therefore, we 
define the small amplitude local wave speed ci oc (x) as 



Cloc(x) = 



1 5s [u(x)] 
8{d x u) 



(24) 



This quantity is plotted in Fig. [5]d for the three models, 
demonstrating the appearance of a dynamical lengthscale 
associated with the nonlinear elastic zone in the crack tip 
region, a lengthscale that is absent in the linear model. 
It is a dynamical lengthscale in the sense that it emerges 
as a result of the dynamics of the crack. One possible im- 
plication of such a lengthscale was discussed in [3, [H[ , 
while additional possibilities should be further investi- 
gated. 

In summary, in this section we have presented a de- 
tailed comparison of the steady state solutions of the lin- 
ear and nonlinear models. We demonstrated that the lin- 
ear approximation breaks down inevitably at sufficiently 
large velocity, if there exists a range of validity for that 
approximation at all. We showed that the linear and 
nonlinear models differ in their limiting velocities, their 
near tip strain, stress and energy distributions and in the 
emergence of a dynamical lengthscale associated with a 



nonlinear elastic zone. In the next section we focus on 
the response of the cracks in these models to small per- 
turbations out of steady state. 



IV. RESPONSE TO PERTURBATIONS: 
LINEAR STABILITY ANALYSIS 

The results presented up now indicate that the lin- 
ear approximation breaks down inevitably at sufficiently 
high velocities. All these results were restricted to steady 
state conditions. However, as mentioned above, one of 
the great theoretical challenges in the field of fracture 
mechanics is the understanding of the origin of crack tip 
instabilities [J, 0|. In this particular respect, the one- 
dimensional model is certainly too simple as we do not 
expect any instabilities to occur here. More specifically, 
since the crack in the one-dimensional model is, by di- 
mensionality alone, restricted to follow a straight path, 
it can at most change its velocity along this predeter- 
mined path. However, it cannot accelerate or decelerate 
significantly due to the global energy balance constraint. 
In contradistinction, the tip instabilities observed exper- 
imentally [l|, [2| involve in an essential way the deviation 
of the crack from the pre-instability straight path. Bear- 
ing this limitation of the one-dimensional models in mind 
and expecting no instability in this framework [Tel . we 
still want to study the response of the steady state cracks 
in these simple models to small perturbations. The mo- 
tivation for that is to gain some insights (or hints) about 
the kind of near crack tip physics that is overlooked by 
the common exclusion of elastic nonlinearities, especially 
as far as perturbations are considered. 

For that aim we perform a linear stability analysis for 
both the linear and nonlinear models. We stress again 
that we do not expect any instability to occur, but rather 
we are interested in the effect of nonlinearities on the 
relaxation time back to the stable steady state. We start 
by defining the coordinate transformation x' = x— x t i P (t) 
and t' = t, where (x,t) is a fixed coordinate system and 
(x',i') is a coordinate system that moves with the crack 
tip x t ip{t). As discussed in [HI EH, sucn a transformation 
is essential in order to avoid irregular behavior at the 
crack tip. Defining u(x'(x,t),t')=u(x,t) and using Eqs. 
^ and ([5]), we obtain for the linear model 

d tt u - 2x t ipd xt u + x 2 lp d xx u - x Up d x u - d xx ii = (25) 
— k 2 (u — 5) — {u + rjdtu — rjxti p d x u) H(l — u) . 

Note that we renamed x' — > x and t' — > t for notational 
simplicity. We are now interested in the time evolution 
of small perturbations of amplitude e around the steady 
state crack tip location 

x U p(t) = -vt - ee ut , u(x, t) = u{x) - eu{x)e ut . (26) 



Substituting these expressions into Eq. 
ing in e, we obtain 



and lineariz- 



(2vuid x + a; 2 ) (d x u — u) — (1 — v 2 )d xx u + k 



2 u 
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— (ii + rjvd x u + Tjcju — r)ujd x u) H{—x) . (27) 

Note that u{x) is simply the steady state solution given 
in Eq. (fTO)) and that Eq. (|2T[) admits a trivial solution 
with u = and cj = 0. This solution corresponds to a 
translation of the steady state solution and is of no in- 
terest here. 

Using u(Q) — 1 (recall that x — is still the crack tip 



location) and u(0) = l, we obtain u(0) — 0. Furthermore, 
substituting u = d x u + u into Eq. (|27[) . we obtain a sim- 
pler problem for u, with u(0) — —d x u(0). The resulting 
problem can be rather easily solved following a similar 
procedure to the one employed in solving for the steady 
states. Specifically, by demanding that w(0) is continu- 
ous, we obtain 



-ArjK 2 [4(1 - v 2 ) + t? 2 (1 + v 2 )] + 2^(4 - r) 2 )^ - v 2 ) [r] 2 v 2 + 4(1 + n 2 )(l - v 2 )} 



4 )2( w 2 



16t] 2 k 2 



(28) 
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FIG. 6: (Color online) Main panel: The relaxation time r 
vs. S /S c for the linear (solid line) and nonlinear softening 
(dashed line) models with k = 0.3 and rj = 0.3. Inset: Zoom in 
on the low S/S c region, where the relaxation time for nonlinear 
stiffening model is added (dashed-dotted line). 



In the range of interest, i.e. for small 77, we find that 
lu < 0, implying linear stability as expected. This result 
implies that the perturbation of the tip location relaxes 
with a typical timescale r — \ui\~ 1 . In the limit /t,?7<C 1 
we obtain r~ (rjK)^ 1 , which is independent of v. 

In order to compare the relaxation time in the linear 
model to relaxation times in the nonlinear models, we 
repeat the linear stability analysis for the latter. Using 
the same notation as before and Eqs. ©-([7]) with 

%tip{t) = ±vt±te u±t , u(x,i) — u(x)±eu(x)e UJ±t , (29) 

we obtain 

=p (^2vuj ± d x ± u) 2 ^ (d x u — u) + k 2 u — 



1 



2 C^QZ ^ *tjb 



2d x ud, 



Z) — ■ 3(1 + cu) 3 (i + d xU y 

u ± r\vd x u + rjuj ± u — r)uj ± d x u) H(^fx) . (30) 



This is the counterpart of Eq. (|2"T|) . Note that the ± 
in Eq. (f2"9"| correspond to the stiffening and softening 
models respectively. Eq. (|30[) is solved numerically using 
a method similar to the one used to obtain the steady 
state solution. The relaxation times are given as before 
by |w ± | . The results are summarized in Fig. [SJ where 
the relaxation times for the linear and nonlinear mod- 
els with k — 0.3 and 77 = 0.3 are plotted as a function of 
S/S c . The results presented in the main panel show that 
the relaxation time in the nonlinear softening model is 
larger than the relaxation time in the linear model. The 
slower dynamics in the nonlinear softening case may be 
attributed to the smaller local wave speeds in the near 
crack tip vicinity. In the inset we focus on the small S/5 C 
where smooth steady states for the nonlinear stiffening 
exist, cf. Fig. [3] The results show that the relaxation 
time in the nonlinear stiffening model is smaller than the 
relaxation time in the linear model. The faster dynam- 
ics in the nonlinear stiffening case may be attributed to 
the larger local wave speeds in the near crack tip vicinity. 
The main result obtained here is that the typical response 
timescale near the crack tip is affected by elastic nonlin- 
earities. This timescale may be of prime importance in 
understandurg the experimentally observed crack tip in- 
stabilities [H, @|- Moreover, this dynamical timescale can 
be interpreted as introducing inertia-like effects into the 
crack tip dynamics [l6| , effects that are missing in LEFM 



V. CONCLUDING REMARKS 

In this paper we investigated the role of elastic nonlin- 
earities in simple one-dimensional models of fracture. We 
were mainly motivated by the recent experimental and 
theoretical findings of [3, H[ that demonstrated explicitly 
the importance of elastic nonlinearities for understand- 
ing the structure of the deformation near a moving crack 
tip. Our results show that the common linear elastic ap- 
proximation breaks down at sufficiently high propagation 
velocities, if there exists a small velocities range of valid- 
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ity at all. This finding is in complete agreement with the 
results of [3, [f| , where it was shown that at high veloc- 
ities the linear elastic approximation of LEFM provides 
un-physical and qualitatively different results compared 
to the nonlinear theory. 

The breakdown of the linear approximation manifests 
itself in marked differences in the propagation velocities, 
including the limiting crack velocity, as well as in the 
stress and strain distributions in the crack tip vicinity. 
The near tip deformation, that is markedly different from 
the linear clastic prediction both in our simple model and 
in Refs. 0, Q , may have a role in determining the stabil- 
ity of the crack tip against perturbations. In this regard, 
we demonstrated the existence of a lengthscale that is as- 
sociated with the nonlinear elastic zone surrounding the 
crack tip. This lengthscale was shown in [5| to coincide 
with the wavelength of the oscillations observed in 
This finding can potentially explain the emergence of a 
non-geometrical lengthscale that is missing in the stan- 



dard approach of LEFM 8] . Furthermore, by studying 
the response of the crack tip to perturbations we showed 
that elastic nonlinearities affect the crack tip local re- 
sponse timescale. This emerging timescale can be inter- 
preted as effectively attributing inertia-like properties to 
the crack tip [IB] , in contradistinction with LEFM where 
the crack tip is regarded as "massless" Q. 

The simple one-dimensional models considered in this 
work offer some insights about the possible importance 
of elastic nonlinearities in answering the long standing 
question of "how things break?" However, in order to 
obtain concrete predictions related to the experimentally 
observed instabilities, one should study crack propaga- 
tion with near tip elastic nonlinearities in higher dimen- 
sions. The novel asymptotic nonlinear solution presented 
in [5| may serve as a promising starting point. 
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